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PROGRAM  TO  CONVERT  NAVIGATION  AND  BATHYMETRIC 
DATA  ON  TAPE  INTO  GEODATA  FORMAT 

1.0  INTRODUCTION 

1.1  Title 

Program  to  convert  navigation  and  bathymetric  data  on  tape  into  GEODATA 
format. 

1.2  Identification  Name 
CONBATH. 

1.3  Classification  Code 
None. 

1.4  NRL  Research  Computation  Center  Identification  Number 
None. 

1.5  Entry  Points 
CONBATH. 

1 .6  Programming  Language 

Language:  3600/3800  Fortran. 

Routine  type:  program. 

Operating  System : Drum  Scope  2.1. 

1.7  Computer  and  Configuration 
CDC  3800. 

1 .8  Contributor  or  Programmer 

Marilyn  L.  Blodgett,  Code  8122MB,  Applied  Ocean  Acoustics  Branch,  Acoustics 
Division,  written  for  the  Environmental  Sciences  Group,  Acoustics  Division. 

1 .9  Contributing  Organization 

NRL  — Naval  Research  Laboratory,  Washington,  DC  20375. 


Manuscript  submitted  December  27, 1978. 
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1.10  Program  Availability 

If  supplied  with  a magnetic  tape,  the  Environmental  Sciences  Group,  Acoustics 
Division,  will  make  a copy  of  this  program  available. 

1.11  Verification 

This  program  has  been  used  and  tested  by  the  Environmental  Sciences  Group, 
Acoustics  Division. 

1.12  Date 
April  1978. 


2.0  PURPOSE 

2.1  Description  of  the  Routine 

The  program  reads  a data  tape  from  an  oceanographic  or  geophysical  cruise  and 
converts  it  into  the  GEODATA  format.  This  is  the  format  recommenaed  by  the 
National  Research  Council  of  the  National  Academy  of  Sciences  with  one  slight 
modification  for  the  navigational  data.  There  is  one  logical  record  (of  80  charac- 
ters) for  each  data  point. 

For  navigation  the  program  converts  the  Julian  day  to  month  and  day  and  con- 
verts the  latitude  and  longitude  degree-and-minute  values  into  their  decimal 
equivalents. 

For  bathymetry  the  Julian  day  is  converted  into  month  and  day,  the  minutes 
and  seconds  are  combined  into  their  decimal  equivalents,  depth  in  uncorrected 
meters  is  converted  to  uncorrected  fathoms,  and  corrected  depth  in  meters  is 
calculated  using  the  Matthews’  table  coefficients. 

The  converted  data  for  both  types  of  data  is  written  on  to  a new  tape  and/or 
printed  out  on  the  standard  printer  (logical  unit  61). 

2.2  Problem  Background 

A program  was  needed  to  convert  the  navigation  and  bathymetric  data  into 
GEODATA  formatted  tapes,  so  that  the  data  collected  could  be  easily  ex- 
changed with  other  facilities  and  so  that  the  scientist  could  use  several  other 
programs  using  GEODATA  tapes  for  input. 
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3.0  USAGE 

3.1  Calling  Sequence  or  Operational  Procedure 
Not  applicable. 

3.2  Arguments,  Parameters,  and/or  Initial  Conditions 
Not  applicable. 

3.3  Space  Required  (Decimal  and  Octal) 

3.3.1  Unique  Storage: 

2740  octal  (1504  decimal)  locations  exclusive  of  system  library  functions. 

3.3.2  Common  Blocks: 

None. 

3.3.3  Temporary  Storage: 

None. 

3.4  Messages  and  Instructions  to  the  Operator 
None. 

3.5  Error  Returns,  Messages,  and  Codes 
None. 

3.6  Informative  Messages  to  the  User 
None. 

3.7  Input 

The  program  has  a Navigation  Parameter  card  which  supplies  the  information 
needed  for  the  GEODATA  format,  the  number  of  files  of  navigational  data  and 
the  type  cf  output  desired.  The  Bathymetry  Parameter  card  serves  the  same  pur- 
pose for  the  bathymetric  data.  This  is  followed  by  a Matthews’  Zone  Values 
card  (or  cards). 

The  navigation  and  bathymetric  data  is  read  in  via  magnetic  tape.  The  formats 
for  both  types  of  data  are  given  in  Appendix  A . Appendix  B is  a complete 
description  of  the  input  setup. 
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3.8  Output 

The  program  will  write  a new  tape  in  GEODATA  format.  Appendix  C shows  the 
GEODATA  format  for  both  navigation  and  bathymetry.  There  is  also  an  option 
for  listing  the  records  on  the  standard  printer  (logical  unit  61).  Appendix  D 
presents  sample  output  listings. 

3.9  Formats 

Appendix  B,  which  shows  the  program  deck  structure,  describes  the  formats. 

3.10  External  Routines  and  Symbols 
ABSF,  ENDFILE,  MOD. 

3.11  Timing 

The  time  required  depends  on  the  length  of  the  input  tape. 

3.12  Accuracy 

Not  applicable. 

3.13  Caution  to  User 
None. 

3.14  Program  Deck  Structure 

Appendix  B describes  the  program  deck  structure. 

3.15  References  — Literature 

M.L.  Blodgett  and  J.V.  Massingill,  “A  Program  for  Storing  Oceanographic  Data 
on  Magnetic  Tape,”  NRL  Report  7861,  Mar.  1975. 

M.L.  Blodgett  and  J.V.  Massingill,  “Program  to  Plot  an  Annotated  Track  or  a 
Track  and  Bathymetry  or  Magnetic  Profile  on  a Mercator  Projection,”  NRL 
Report  7930,  Feb.  1976. 

L.  LaLumiere,  “Program  OCEANO”  (not  published). 

4.0  METHOD  OR  ALGORITHM 
Not  applicable. 
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5.0  FLOW  CHART  AND/OR  SOURCE  LANGUAGE  LISTING 

The  flow  chart  is  given  in  Appendix  E,  and  the  listing  is  given  in  Appendix  F. 

6.0  COMPARISON 

There  are  no  other  known  programs  available  for  comparison. 

7.0  TEST  METHODS  AND  RESULTS 

Samples  of  the  GEODATA  format  for  both  navigation  and  bathymetry  are  given  in 
Appendix  C.  Sample  output  listings  are  given  in  Appendix  D. 

8.0  REMARKS 

Program  CONBATH,  in  order  to  convert  the  bathymetric  data  into  the  GEODATA 
format,  requires  the  different  Matthews’  zones  the  ship  passed  through.  One  method 
to  obtain  this  information  is  as  follows: 

• Run  Program  CONVNAV  to  convert  the  navigational  data  into  GEODATA 
format. 

• Using  the  GEODATA  formatted  tape,  plot  the  navigational  data  with  every 
fifth  fix  annotated  on  a Mercator  projection. 

— Use  Program  OCEANO  (see  references)  with  a scale  of  0.3966. 

— Use  Program  MERCATOR  (NRL  Report  7930,  see  references)  if  the 
required  heights  can  be  measured  from  an  already  existing  chart. 
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Appendix  B 
DECK  ASSEMBLY 


3 Program  CON  BATH  (binary  deck).  (If  the  Fortran  source  deck  is  used  instead 
of  the  binary  deck,  a Fortran  card  is  required  after  the  Equip  cards.  In  addition, 
a Scope  card  and  a Load  card  must  follow  the  source  deck.) 


2 Equip  cards  (logical  units  15  and  40) 


1 Job  card 


Card  Column 

Number  Title  Number 


Description 


1 Job 


2 Equip 


3 Program 


1-21  7/9  JOB,  Charge  No.,  ID  No.,  time.  See  page  2-2  of 

the  3600/3800  Computer  System  Drum  Scope 
Manual. 

1-18  7/9  EQUIP,  15  = **,  RO,  HI 

7/9  EQUIP,  40  = **,  WO,  HI 
15,  40  = logical  unit  number; 

RO  = read  only; 

WO=  write  only; 

HI  = high  density. 

Deck  of  This  is  the  main  program  and  associated  subroutines 
CONBATH  in  binary  form. 


4 Run  1-13  7/9  RUN,  T,  P,  R,  M,  D 

T *-  time  limit  in  minutes; 

P = maximum  number  of  print  or  write  operations; 

R,  M,  and  D may  be  left  blank.  See  page  2-15  of  the 
3600/3800  Computer  System  Drum  Scope  Manual. 


BLODGETT 


Number 

Card 

Title 

Column 

Number 

Description 

5 

Navigation 

Parameter 

1+2 

75  (IYEAR) 

Last  two  digits  of  the  year. 

3-10 

A64-17-8  (CRUISE) 

Cruise  identification.  Any  eight-digit  alphanumeric 
number  can  be  used. 

11-15 

1 (NOFILE) 

Number  of  files  of  navigation  data  on  the  input  tape. 

16-20 

0 or  1 (IOUT) 

0 = just  write  the  navigation  data  in  GEOD  .TA 
format  on  the  tape; 

1 = write  the  navigation  data  in  GEODATA  format 
on  the  tape  and  also  print  it  out. 

6 

Bathymetry 

Parameter 

1+2 

75  (IYR) 

Last  two  digits  of  the  year. 

3-10 

A64-17-8  (CRUISE) 

Cruise  identification.  Any  eight-digit  alphanumeric 
number  can  be  used. 

11-15 

1 (NOFILE) 

Number  of  files  of  bathymetric  data  on  the  input 
tape. 

16-20 

8 (NOMZ) 

Number  of  Matthews’  zone  changes  during  the  ship 
cruise. 

21-25 

0 or  1 (IOUT) 

0 = just  write  the  bathymetric  data  in  GEODATA 
format  on  the  tape; 

1 = write  the  bathymetric  data  in  GEODATA  format 
on  the  tape  and  also  print  it  out. 

7 

Matthews’ 

Zone 

Values 

2-5 

0705  (ITIME  (1)) 

Time  before  or  when  the  first  data  point  was  taken. 
The  first  two  digits  (CC  2+3)  specify  the  hour  and 
the  last  two  digits  (CC  4+5)  specify  the  minutes. 

8-10 

262  (JUL  (1)) 

Julian  day  when  the  first  data  point  was  taken. 


8 


NRL  REPORT  8306 


Card 

Number  Title 


Column 

Number  Description 

14+15  2 (MZ(1)) 

Matthews’  zone  the  ship  was  in  when  the  first  data 
point  was  taken. 

17-20  1215  (ITIME  (2)) 

Time  when  the  ship  passed  into  a different  Matthews’ 
zone. 

23-25  262  (JUL  (2)) 

Julian  day  when  the  ship  passed  into  a different 
Matthews’  zone. 

29+30  3 (MZ  (2)) 

New  Matthews’  zone. 

32-35  1956  (ITIME  (3)) 

Time  when  the  ship  passed  into  a different  Matthews’ 
zone. 

38-40  263  (JUL  (3)) 

Julian  day  when  the  ship  passed  into  this  Matthews’ 
zone. 

44+45  4 (MZ  (3)) 

New  Matthews’  zone. 

47-50  1000 

T.  me  when  the  ship  passed  into  a different  Matthews’ 
zone. 

53-55  264 

Julian  day  when  the  ship  passed  into  this  Matthews’ 
zone. 

59+60  3 

New  Matthews’  zone. 

62-65  1300 

Time  when  the  ship  passed  into  a different  Matthews’ 
zone. 

68-70  264 

Julian  day  when  the  ship  passed  into  this  Matthews’ 
zone. 
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Card  Column 

Number  Title  Number  Description 

74+75  2 

New  Matthews’  zone. 

As  many  Matthews’  Zone  Value  cards  as  required 
may  be  used.  The  number  of  change  sets  must  equal 
columns  16-20  of  the  Bathymetry  Parameter  card. 
The  last  time,  the  Julian  day  and  Matthews’  zone 
stated  must  be  on  or  after  the  time  when  the  last 
datum  was  taken. 

8 End-of-file  Terminates  the  run. 


Appendix  D 

SAMPLE  OUTPUT  LISTING 

I 
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Appendix  F 

SOURCE  LANGUAGE  LISTING 


PR0G«AM  C0NVF.  ATM 
C£0MEN/l(-)F.RH/lG6F 

C I A THfc  feVtNT  EF  AN  l/E  CGNVtFSIEN  ERK£R  AT  STATpMf^T  209  Ctlfc  T p 

C A A INVALID  CHARACTER,  FUCGF  WILL  SFT  in£F  T 3 l 

CALL  fudge 
lE^O 
IF  lx  = 0 
ITM2Ae*c 
UHil 
I CNT  *o 

RfcAD(60,200)!YR,ChlISE,NEFTLE  . IELT 
20C  FERMATtl2,Afl.2!5> 

3CC  READ(1‘>.100)lLI:V,FR.XKIN,XI  Ar,XLAfSXL6D,VL0H.  ITYPE.  I f I X 
IOC  FERMAT!  l3,F2,F4.i,F3,f5.2,F4,F5.2, 12, J4) 

IF  < IflCwfcCK,l5)  30C.3C1 
331  IF<E0F»15)  3C2.303 

C CONVERT  JULIAN  CAY  T£  MEATH  ANE  CAT 
303  IF  < IG£F  ,Ea,l)  GS  Tf  444 

CALL  JULtAN(ITR,JLIT,IC,IM,LFTR> 

IF  ( [F  [X.Lfc,  IFU2)  GE  TE  2 0C 
IF IX2«IFIX 

xmin*xmin»io 

C CHANGE  LATITUDE  ant  LENGJTLDE  TE  DECIMAL 

MiNus»40Joooccoooocccoe 

XLAM  = ARSF(Xlam  ) 

XL0M*AHSF(XlFM) 

KEYeXlAD.aND.mmjs 
IF  (KBr.tO,MIMS)27,2f 
27  CEMlAT*XLAD-Xlam/4C  , C 
GE  TP  29 

2f>  C£MlAT»XL*D  ♦ XLAM/6C.C 
29  kETeXlRD.  AND. MNUS 

If  <KtY.fcO.MINLS)30,3l 
3C  C£MLPN*XLen-XLf.M/tC.C 
GE  TE  4C0 

31  C£MLGN«XL9n  * XL^M/bt.C 

40C  MMTE<4C.40l)CRLlEF.!TM2NE,iTR.IM,ID,HR,XMlN.CPMLAT,CdMLPN,ITYP6, 

1 1 F I X 

4 C 1 FERMAT <A6, IS. 312, 1XF2 ,F3,F » .4,FS ,4, 12,7X15. 24X> 

IF <IEuT.Nfc,l)  GE  IF  3 0C 
ZMIN«XMIN/10,C 

IF  <MED<Niim,60>.NE.1  ) GE  TS  6 CC 

501  lR I Tt (61 , 501 ) 

5C1  F£HMAT<1H1,123MSH1F  ANt  CRlIEE  TIMF  yeah  MPNTh  DAY 

1 H«UR  MIMLTE  LATITl.CE  LSNu I TUQF  FIX 

2 FIX) 

NR 1 TE (*1,502  * 

502  FEHMATdH  ,24HC&ATJFICAT1EN  20Nfc,7OX,22HDESCRIPTieN  NUMB 

1£R ) 

NLM*1 

60C  V.fiITE-(61,116)CRLISF,lTM2Nb.  IYR.  JM,  lo.HH.ZMl  Y.CPMLAT.CBMLPN,  ITYPE, 
1IFIX 
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■ GE  IS  BEST  Q.VUil'** 

. t°DUC 


£ I U*  V - 


BLODGETT 


116 

F6RMAT(1H  ,2*AP,7xI5, 
1,10X15) 

7>I?,7XI2.7XI2,7XF2,7XF4,1,6XF8,4,5»FQ.4,<)XI2 

NLM*NUm  ♦ 1 

GE  TB  300 

444 

1 6EF  s 0 

GE  T£  300 

302 

ICNT=1CNT  ♦ l 

IF < ICNT.EQ.NEF IlF  ) GE 
GE  T0  300 

T 0 7C0 

70C 

ENDFILF  40 

CALL  BATHY( ITMZNE) 
ST0P 

END 

Itsnt  cpnvbath 

PRCGRAM  LENGTH  0 C A A 7 

ENTRY  PEIMTS  C 9 N V P A T ► 0C135 
BlPCK  NAMES 

I HERR  CC0C1 

FXTEPNAl  SYMMOLS 

OdCtNTRY 

ThFNP, 

cbcstgps 

OoUlICT. 

FutGF. 

Julian 

ratpy 

xmecf 

ObCIFEMF 

Ob^lFIRC 

EFT, 

TSF, 

STM. 

fJNFlNGL, 

00125  SYMBOLS 


SL3R6UT 1 Nfc  jLLIAM  f I V » wOt V . I C 

.IM.LFYR) 

F 

1 

ENVERTS  JULIAN  DATF  I NTG  D* Y , 

VEN I K , YEAR 

F 3P 

ANM0TAT1PN  IN  PCFAN0 

F 

?' 

IY  s YEAR,  JLtY  * „LL  1*N  C»Y, 

IG  * CALFNrtR  PAY,  I M B CAlEMDER  MSN 

F 

3 

EMPlUiENTS  0F  P0B  FECbN  - CECE  fcl74  • 17 

NE>V 

1971 

F 

4 

M £ D I F I F.  D BY  LEEN  LA  LLMlFRF  • 

reck  0i7"  - 

17 

N0V  1971 

F 

r 

5 

A 

riMENSISN  J J ( 1 3 ) , A A ( 12  ) 

F 

7 

TYPE  JMTEGFH  AA 

F 

* 

data  ( ( J J ( I ),  I - 1,  13)  S 0, 

31,  59,  «0. 

120 

, 151,  1«1,  212.  743. 

F 

9 

1273,  304,  334,  3*5) 

F 

10 

CATA  ((AA(I  >.  I s 1,  12)  * 1, 

2 , 3 , 4 , > , 

6, 

7,  ",  9.  10.  11.  12) 

F 

11 

LL*  I Y 

F 

17 

MM  * JLD Y 

F 

13 

IF  (LL.LT.O)  GE  T £ 1C 

F 

14 

IF  (LL.EO.O)  GE  T£  1 

F 

15 

IF  (M(ji)(LL,4).EC.0)  CE  TE  4 

F 

1* 

1 IF  (Mm.LF.0,FR.MM,GT.3<5)  G6 

TE  4 

F 

17 

2 CE  3 K«2 , 13 

F 

1" 

IF  (Mm.LE, JJ'M  ) ) GE  TE  6 

F 

19 

3 CENTIME 

F 

20 

18 


JJUS  PAGE  IS  BSSX  HUALin 

Qyr'l  /UWUSHBD  TO  DDC  - 

NRL  REPORT  8306 


4 LFYr*2 

IF  (MM.LE.O.RR.MM.CT.SOPI  GO  1 ? 10 
IF  (MM.LE,5v>  GO  TP  ( 

IF  (MM.EO,60>  Gt)  TP  5 

HNSHM-I 

go  tg  2 

5 11*2? 

K*3 

Gt  TP  7 

6 I I*MM-JJ(K-1) 

IF  (LL.EO.OJ  OP  Te  6 
IF  (M0’l(LLi4>.Nfe.O>  CO  TP  f 
IF  (P0t)(LL.4>,FC.O.APC,MK.LE,5<!)  GO  TO  8 
KMMK*1 

7 IF  (LL.GE.10>  Ge  TP  F 

8 I C * I I 

IhsAA(K-l) 

9 RfcTuHN 

10  PRINT  11,  LL.8>' 

GO  T«  9 

11  FORMAT  (lX**YfcAR  s • , 12 # 1 X . '.LI  1 AN  DAY  s * , I 3 , 1 X . • A RF  RAD  DATES*  > 
FND 


F 21 
F 2? 

F 73 
F 2* 

F ?S 
F ?ft 
F 27 
F 28 
F 29 
F 30 
F 31 
F 3? 

F 33 
F 3A 
F 38 
F 36 
F 37 
F 38 
F 39 
F 40 
F 41 
F 42 
F 43- 


5.4DS  LILIAN 


ITEM  JULIAN 

PROGRAM  LENGth  3C3C1 

ENTRY  POINTS  JULIAN  0CC83 

EXTERNAL  SVMU0LS 

thrnd. 

08LLICT. 

XMPGF 

STM, 

0N51NGL. 

0C116  SYM"0LS 


I C E N 7 FUDGE 

PfifiGRAM  LENGTH  GC073 

entry  points  FurGE  ocooo 
block  names 


I0ERR  OCOOl 

external  symbols 


ELC. 

06CERSET 

EXT 

ELD. .080EPSFT 

ENTRY 

FUDGE 

00000 

I0ERR 

bleck 

oooco 

COMMON 

lBfcF 

ocooo 

00 

G 

OCCUO  FLCGF 

6 C 7 

0 

00 

0 

00C00 

0C0C1 

77 

2 

OCCOO 

csta 

xn 

20 

0 

P00C15 

00002 

86 

1 

POCC11 

S 1 L 

RFST&H.l 

80 

0 

OCCOO 

19 


1S  BESX  QUALITY 

i?o!«  rswns^T0DDC  " 
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00003 

77 

1 

0*0  30 

fc  A A 

OFUUGF.l 

10 

0 

P0CC17 

00004 

50 

1 

OCCOO 

EAI 

0.1 

50 

0 

oocoo 

0C0C5 

6 3 

•> 

10017 

SFYT, 

A0.E15  x,.l 

50 

5 

P00C14 

0C006 

50 

1 

0CC17 

E A I 

15,1 

12 

3 

POO  C 02 

LCA 

F i DRE*2 

00007 

*3 

0 

14C03 

iFYT  , 

A3?,E3  X.,1 

50 

b 

P0CC14 

CC010 

12 

0 

PO0C14 

LCA 

X 

20 

o 

X77/77 

STA 

C RqpRSE  t -? 

CC011 

50 

1 

77777 

RESTtlR 

fc  A I 

•*.i 

50 

0 

oocoo 

0001? 

77 

2 

OCCOO 

CLCA 

x*i 

12 

3 

POCClS 

CC013 

75 

n 

POOCOO 

si* 

fudge 

50 

5 

OCCOO 

00014 

63 

c 

OCCOO 

X 

LP*F 

• • 

00 

j 

77777 

C 0 015 

00 

0 

OCCOO 

£ C T 

0 

DC 

3 

050C0 

0001* 

00 

3 

OCCOO 

ect 

3 

00 

C 

oocoo 

00017 

77 

2 

OCC  JO 

FLEGE  ,1 

ecu 

x*l 

20 

} 

P00C15 

00020 

10 

j 

OOCOl 

EAA 

1 

20 

n 

COOCUO 

STA 

IPEF 

C 0 021 

77 

2 

OCCOO 

CLCA 

X*1 

12 

O 

P0CC15 

0002? 

63 

3 

OCCOO 

LE*p 

(5)FLU,*10 

01 

r 

J 

X77777 

fcAC 


00007  S^OLS 


A 


SLPRPliTINfc  B*Tt-Y(  I TK7NF  ) 

ClMENSlbN  lTIf'F(40),*LL(40>,»'Z(40) 

KCM«1 

TCNT*0 

ITM2Nfc«0 

RfcAD(60,100>IYR,CRLlSt,NEFUfe.A£»'2,IGl,T 

ICC  F£Rmat( I2,*b.3I5> 

RfcAD<60,40C)(IT[ME<l)t»liUl)»l'7<l>.l«l»NPHZ) 
4CG  FeRM*T(l5I5,5x> 

*AT«PZ«1> 

!*? 

300  RE *0(15,200)  «trY.>R. AKIN, SFC.AL *T .XlON.UNCOR 
200  F£RhAT(  J3»3F?,?f<,4,F7.0> 

IF  ( Jt!C*fcCK.l5  )300,30l 
3C1  IF (EPF, 15)302.303 
C CflKVERT  JULIAN  C4V  Te  KfATH  *» t t*T 
303  if  (jcef.eo.i)  r.e  tf  <4* 

CALL  JULlAN( IYR, JLTT, It  , IM.LFYR  ) 

C CHAAGfc  A i NyTfcS  AAC  SfcCEAl*  Te  EECIPAL  ECOIVALtYl 

»<iNus«4c:ooootooooccccF 

SE  C*  AgSF (StCI 
K£Y«*MlN. ANU.R IAUS 
IF (H£Y.t3,NlNLS)27,2E 
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27  CeMLAT»XMIN'-Str/60  . o 

ge  re  ?9 

2®  r£MLAT*XMiN  * SEC/zG.O 
C CAtCILATfc  CORRECTED  KETFKS 
2?  C6pTh*UNCRH/1 ,82fcf 
ICL0Ck«hR*1O-C  . XmIN 
ICEPThsDEPTh 

IF  (JUDY.FQ,  JLL(  l ) ,ANC  . ICLEiCK  ,GF  ,lTtMt<l>)  3«  T0  600 
GE  TF  AC2 
60C  KATiMZU) 

1*1  ♦ 1 

632  CALL  RTCPRI  rEPTH.FAT ,F6CFaT ,wFRt ET.METuNC, MTQC > 

C6MlAT»GHMLaT«io 
vtfcP7H«UEpTH  • 10 

WfilTfe<4C,605KRLlSEiITW*KE,  lYfi.lF  , I D . ► » , C0Ml  A T , XL  A T , XLON  , M"EP  TR , 

1KGRME  T i MAT 

6G5  FERMAT < AS, 15. J 12, 1>FJ ,F3,FF ,4,rS ,4,10X215, 12, l6x> 

IF  (ieuT.NE.1)  G0  T £ 300 
ZMN  = C0MUAT/1C.O 
IF(m£U(NUM,6C).NE.1  ) (if  T 0 7C0 
M Mi  1 

wp:tf(6i,5c«) 

WRITE (61.5C5) 

5C4  F£Rmat<1Hi,i33MSMIF  A NC  CRLlSE  TIME  YEAR  MRNTh  CaY  hPUR 
1 M^UtE  LATITUCE  L0F  Cl  tLEE  LnCG'RRFCTED  UNCSRRECTEp  r0«fl 
2FCT£P  MATTHEWS) 

505  F6RMATI1H  ,22HDEUIFICATIFn  ZeN£,65x.44WFATncmS  >'ETERS 

1 METERS  ?GaE> 

70 C WRITE (61,70l)CRLlSE, ITM2AE. IYR. IF , ID,HP,7M1R.XL*T,Xl0N,  PERTH, 
1LNC0R,K0RmFT,FAT 

701  F£RmAT(1h  ,2XAfl,5xl5,5*I?,5XI2.5xi2,5xF2,5XF4,l,4XF8.4,3xF9.4, 

1 5XFfc,l,SXF6, 1,7X15,9X12) 

NLMiKUM  ♦ 1 
G£  T6  300 

3C2  ICNT*ICNT  ♦ 1 

IF ( ICnt  ,EO,N£f  ILE  ) G£  T0  bCO 
G£  T£  300 
444  I £EF  *0 

G£  T£  30C 
eoc  erofile  40 

REWjRi)  4'.' 

REWIND  15 
RET ijRn 
EM) 


KENT  fcATMY 

PROGRAM  LEAGTH  0C753 

ENTRY  PElN^s  Matey  0C334 

EXTERNAL  SYMt<0LS 

thfnc. 

oir.ioioc 

ohgcict. 

JULIAN 

MTCES 

xmel'f 

OttGIFE^F 

'jaciFjoc 


nnnnnnnonnannnnonnnoonnnnnnn 


BLODGETT 


C 


SLBREyT ] NE  MTCSR  ( ICEP.MT.KCCF  .KCEH.METuNC.iTDC)  J 

J 

JAN  10.  19^2  J 

J 

H6D  JF  JED  FSR  CCC  3800  EV  LFEN  lA  LLmIERE  • CODE  817*  • 20  MAR  1972  J 

J 

SLBRBUT  I N£  MTCeRCAFG>,  MATTHEWS  TABlE  CORRECT  I BN  J 

VFRSJ0N  WHERE  C8EFFJCIENTS  ENTERED  AT  RUN  TIME  J 

AtCULATES  cbrrectec  CEFTW  IN  PETERS  WHEN  GIVEN  J 

MATTHEWS  TABLE  C6EFFICIENTS  AND  UNC0RRECTED  J 

DEPTH  IN  FAThEMS,  j 

CEPIED  FRBM  WWEl  PR0C-RAH  FER  IBP  1130  J 

J 

INPUT,  J 

IDEP  ■ UNCBRRECTEr  CEPTW  IN  FATH0MS  J 

mt  « Matthews  Table  ceefficientu  - 52)  j 

j 

CLTPUT.  j 

KCDF  « CERRECTEr  DEPTH  IN  F ATt-@MS  J 

KCDM  « C6RRECTEC  CEPTH  IN  PETERS  J 

MTDC  ■ PATHEwS  TAELE  DEPTH  C6RRECT 1 9N  (METERS)  J 

J 

INTERNAL  j 

NUDF  WATER  DEPTH  UNCERRECTED  FATH0MS  J 

WUDF  WATER  DEPTH  UNCERRECTED  F ATH0HS  J 

WUDM  WATER  DEPTH  UNCERRECTED  mETFRS  J 

WC0R  C6RRECTI0N  IN  MFTERS  J 

J 

SETTING  MATTHEWS  TABLE  CCEFF1CIENTS  SUPPLIED  BY  C.  G*NTAR  J 

Cl  MENS l BN  MTN(52),AHT(52),BPT(52>,CmT<52),DHT(52),EmT(52),FMT(52)  J 

J 

REAL  IDEP 

DATA  <mtn«1,2,3, 4, 5,6,7,8,9,10,11. 12,13,14, 15,16,17,18,19,20,21, 22  J 
1,23, 24, 2^,26, 27, 28, 29, 30, 31, 32,33,34,35, 36, 37,30,39, 40, 41,42,43, 44  J 
2,45,46,47,40.49,50,51,52)  J 

DATA  (AMT«-,0S, ,2l, .19, ,25, ,35, ,43, ,12, -.21, . 76 , - , 19. * . 07. 1 . 29  , .34  J 
1 » * , 0* , ,12, ,39, 1,18, 1,25. *,11, ,64, «,3l, .43, *, 29, -.0*,*, HP, *,06,*, 2 J 
23, . 0* , 1 , 14 , • , 38 , • 2 , - , 5 6 , ,34, ,58,1,15, ,50, ,01, ,59,2,49, *97,3,93,1,9  J 
36, -.15. ,09, ,06, .42, -.05, . 03 , - , 41 . - , l9, - , gl , ,07)  J 


CAT a (HMT.-, 01603,-. 00247, . 006e85 ,, 000*3?, ,0157i9, .004269, ,013881.  J 
1.025565,  .030789, ,02  3035, .0316  73, .036491, .031972, .0*9780, ,032253. .0  J 
2*0042, .037533, ,0*3541, ,0*6531, .034097, . 019659, , 025093, . 012063, , 005  J 
3833, -,003139, -.006374,-, 00994, -,012 026. -.019317, -.01 1231,, 0341 03, . J 
* 037i2l, ,037074, ,03524, .031429, .037872, .0*3038, ,0395, ,025775, ,03583  J 

50. . 021001. .0186. .021906. .013118. .001109. .0*4699, .02651, ,020006, .03  J 

65524..  03841. ,C*90O7,.OCOe32)  J 

CAT a (CmT.,130695,- .C28055,-. 033*16, ,030785, -.169795. ,049419, -.059  J 

193 1,-, 156579 r-,141286, ,00879 ,- .101806, -.121871,-. 0*0565,-, 233906.  . J 
2011260 ,-. 363995 ,-.227838,-. 239633.., 303032 , '.2*561, -.0*4004 ,-.2165  J 

309.. .014753...0092C9,, 021 74  6, ,062109, ,101101.  .11045,  .193082, .06327  J 

48, -.255508, -.174712,-, 196  492, -.149662, -.141*22, -.17812, -.210*77,-.  J 
5103*53, -.092957, -,142 494, -I  0*741, -,0227i,-, 117*73, -,03’335, ,0236*8  J 
6, -. 271705, -.0* 1812, ,0*8156, -, 059269, -,097981. ,001276. .161859)  J 

CATa  < DmT«-,C* 0587. .£ 17 161, -,00*71, .012525,, 115271, -.01 5997, .0*609  J 
14,  .0*3656, ,05727,-. 01 941, .036137, . 039775, , 00*6*0 .. 083407 030352.  J 

2.207727. . 102774. .095135. .107189. .123539. .022*5, ,119*29. ,000118, .02  J 


1 

2 

3 

4 

5 

6 
7 
0 
9 

10 

11 

12 

13 

1* 

15 

16 
17 

1* 

19 

20 
21 
22 


/ 


23 

2* 

25 

26 
27 
20 

29 

30 

31 


32 

33 

3* 

35 

36 

37 
3* 

39 

40 

41 

42 

43 

44 

45 

46 

47 
40 

49 

50 

51 

52 

53 

54 

55 
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C 

C 


C 


3 42 32. . 01811, -.013313, -,020*12,-. 0305*5. -.0*1797, -.001826, ,121471,. 
4 06  4 336. ,381851, .056314, ,0  55151. ,067, , 083684, . 0201  31  . , 0 3766 1 , . 0 9235 

59. .  01 88, 48,, 0CS24..C6C 13., 0237  95  , ,00 74 13, .103092, .1*79  05, .004  016,  .0 

670316..  09786,, -.00C464.-.162C7P) 

CAT  A tf.MT  1,09  2 39  4,  . 0 4 89  V5  , 127 1 36  , - , 02266  , - , ?49X66  , . 06975  087349 
1. -. 154994 C75775, ,C81 004, -,03567b, -,040077. ,023167. -.115908. .107 
256,-, 47 C13’,-. 16656 8, -,137?7 8, -,366 ob7... 235691, -.029841, -.232583. 
3, 015918, -,04ft 104. -.045777, .024586, ,o37 058, ,062662 ,,110085,. 001*86, 
4-. 21*992, -,0Pefi5,-. 1381 37, -.C79156,-. 07518, -.09377, -.136212,  .01265 
52, -,  041351, -.186534, -,C1*159, , 0 0 1 108  .-. 09774 0256&2 .-, 0 0 *7l7  1 
642 74  8, -.983392, . 00 C 155 , - . 17  0 839 , - . 25?7 . , 74 145 ,, 983392  ) 

CAT*  ( FmT 0 721 7i,-. 15 61 4e,-.  23 102,  .012056,  .190978, -.076 755,  ,0646 

179. . 109779,  ,038823  ,-.0753  41  , , 0 11294 , . 015787 ,-. 027333  , , 0674  95 , -. 096 

2216. . 403122. .10073. .074903. .267798. .169277. .021495. .16507, -.0289 03 

3. . 035913. . 040884, -.01 41 77 ,-,022102,-, 045831, -.06729, .001066, .14100 

47. . 0,968, ,094  7 43. ,0*8  08  5, , 0395 34  ,. 050998 ,. 0B9l76 03255? ,.  01 7 001  , 
5,1*0729, , 00386 7, -.004373, .058355, .010515, .000162, .067775,2.003705, 
6-, 00155, ,1*9654, .2 35294, -.623288,-2,003203) 

NZErE«0 

MDF«IDEP 
IF  (NUOF)  3,5,1 
1 IF  (52-MT)  4,2,2 


J 56 
J 57 
J 58 
J 59 
J 60 
J 61 
J 6? 
J 63 
J 64 
J 65 
J 66 
J 67 
J 6« 
J 69 
J 70 
J 71 
J 7? 
J 73 
J 7* 
J 75 
J 76 
J 77 
J 78 
J 79 
J 80 


2  WlDFMDfcP 

WIDm*WUDF»1,82P8 

WC9R«AMT<Hl)*Ef*T(pT  >»WLD8*fPT  <PT )*1 ,R-04»(^UDM4*2>*DMT(KT)«l.E-07* 

1 ( kUDM*«3  > *EMT  («T  >•  1, E-ll«<  HOP  ••  * >«FPT(MT)»  1.6-15* <VIUDM*»5> 

RCDM»MllDM*WCeR 
KCDm*rCDm*0 ,5 
PTDCewC0R*O,5 
KCDFs(WCDm*0,54681)*C,5 
PfcTgNC*WUDM 

return 

c 

3 PRINT  6 

GE  TO  5 

c 

4 PRINT  7 

C ERR0R  EXIT 

5 KCDM=NZfcR0 
8TDC=NZE80 
KCDF«NZER0 
RETURN 

c 

6 FERmAT  ( IX • *C£PTH  fiEAC  IS  NEGATIVE*) 

7 FERMAT  (1x,*RATTHEWS  ZENE  M.P6FR  IS  GREATER  Than  52*) 

END 


J 82 
J 83 
J 8* 
J 85 
J 86 
J 87 
J 88 
J 89 
J 90 
J 91 
J 9? 
J 93 
J 9* 
J 95 
J 96 
J 97 
J 98 
J 99 
J 100 
J 101 
J 102 
J 103 
J 10* 


ICENT  PTC SR 

PRGGRAR  LENGTH  31015 

ENTRY  PEInTS  MTCER  0C577 

EXTERNAL  STMH0LS 

O1C10100 
ThFNC  , 

38CCICT, 

STR, 

C0C74  SvMii0L5 

LftAD 
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